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ABSTRACT 



Problem areas in the field of computer graphics, as applied to three- 
dimensional space, are introduced through a discussion of hidden-line 
elimination and perspective views. An adaptation of a simple, fast 
algorithm for hidden-line elimination is presented. 

Graphics terminal display capability for a two-dimensional repres- 
entation of a surface is made available. Interactive extensions to the 
basic program are developed to enhance potential applications for 
design and analysis. 
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I. INTRODUCTION 



Man-machine communication of three-dimensional information pertinent 
to shape descriptions has been, in recent years, a subject under investi- 
gation. New techniques to facilitate design and analysis activities have 
been adopted using graphical methods. 

A. THE GENERAL PROBLEM 

The implementation of translation and rotation capability to include 
hidden-line removal for three-dimensional objects in perspective views 
is indeed a formidable problem. The goal, a procedure generalized to 
include a large class of arbitrary surfaces, objects, and space-curve 
intersections, has long eluded investigators in this field. 

The difficulties encountered in this area have been compounded by 
the processing time requirements of generalized procedures. Approaches 
taken in the solution of the general problem have primarily been to 
narrow or isolate one or more sub-problem areas for analysis. 

B. THE SPECIFIC PROBLEM 

In order to limit the extent of investigative efforts in relation 
to the general problem, a more specific goal was developed. Graphical 
display output of a surface or series of curves to include hidden-line 
elimination was the primary problem area to be investigated. A constraint 
on processing time requirements was added for real-time analysis capa- 
bility. Exploitation of the highly dynamic characteristics of inter- 
active computer graphics was also to be investigated. Operator 
modification of the displayed structure and selection of design and 
analysis options were seen as potential candidates for implementation. 
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With generalized objectives in mind, a potential field for applica- 
tion was sought in order to further specify the nature of the techniques 
to be developed. Signal processing provided a representative area for 
possible application of resulting methods. An example of the three- 
dimensional display structure required was contrived. A series of 
curves, each curve to represent a frequency-amplitude function at a 
point in time would be exhibited. In terms of rectangular coordinates, 
the ordinate might represent frequency, the abscissa depict signal 
amplitude, and time as the depth dimension. An example of a processed 
signal represented by such a display structure is shown in Figure 1; 
with separate curves corresponding to the signal at discrete points in 
time. 

To facilitate analysis of the display, hidden-line elimination was 
deemed essential. Visual determination of relative amplitudes among 
succeeding curves would be difficult, at best, with wire-frame 
representation. 

A real-time display environment was part of the problem specification, 
thereby adding a processing time constraint. Relating this condition to 
the signal processing display example, time delays between drawing of 
the curves should be minimized. Since each curve would represent signal 
sampling at a discrete point in time, visual interpolation between curves 
by a viewer would provide the concept of a surface, or continuum of 
sampl ing. 

Rotation of the displayed structure was desired to present the user 
with optional aspects, or viewing angles.. For example, in Figure 1, the 
signal ridge would not be visible if the display were oriented (rotated) 
with receding curves to the left relative to the initial signal curve. 
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Figure Example of Display Structure. 

Perspective views of a family of curves were not initially considered 
for investigation. It was felt that qualitative comparisons would depend 
on visual reference among similar curves, i.e., not transformed to pro- 
duce perspective depth cues for three-dimensional objects. Offset in 
viewing angle by simulated rotation was anticipated to provide the 
necessary time-step (depth) discrimination. 

Some degree of interaction between an operator and a graphics dis- 
play unit linked to a computer was demanded. The potential of a cathode 
ray tube (CRT) for display would not be fully used without regard to 
man-machine communication links. The express form of interaction was 
not specified. However, enhancement of quantitative analysis and various 
display options were areas for investigation. 
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No mention has yet been made of the data format requirements for 
the display. The intention was to provide a highly useful output 
medium. Any preprocessing of data to provide a suitable format for 
manipulation and graphical display was assumed. 

C. PREVIOUS WORK 

Following are representative summaries of research to date in the 
field of object graphics display. Hidden-line removal and perspective 
view considerations are emphasized. 

Johnson [Ref. 8] presented a report including a good background for 
his then present activity, the manipulation of straight-line, wire-frame 
figures in three-dimensional space. The visual presentation .consisted 
of three orthogonal views in addition to a perspective view of an object. 
Johnson states that the two-dimensional views without perspective would 
not convey the proper depth information. The perspective view alone 
also did not visually display all the correct information, particularly 
with regard to wire-frame representation. Johnson also chose to emphasize 
the importance of real-time, bilateral communication, as dfd Sutherland 
[Ref. 21]. The immediate future problems, as they appeared to Johnson, 
were allowance for arbitrary surfaces, surface intersection determination, 
and hidden -line elimination. 

An article by Puckett [Ref. 20] considered rotation and translation 
of three-dimensional objects prior to perspective projection on a two- 
dimensional surface. Stereographic pair applications were also discussed. 
Again, wire-frame representation imposed limitations on the usefullness 
of the display. 
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More recent work, primarily concerned with perspective representa- 
tions, has come to include consideration of hidden-line removal as an 
integral part. 

Kubert, Szabo, and Giulieri [Ref. 9] introduced their investigations 
of functions of two variables by discussing contour plots and planar 
projections. Contour plots had an inherently quantitative nature. 

Planar projections were used where qualitative descriptions were essential. 
Families of curves used to define an object were subjected to prospective 
transformations.. Criteria for the visibility of points defining an 
object were established. 

Galimberti and Montanari [Ref. 5] had shifted the emphasis to hidden- 
line elimination in their investigations of perspective display of con- 
cave and convex plane-faced objects. Computing time proved to be a 
major constraint. 

An earlier article by Appel [Ref. 1] included a discussion of 
previous research into hidden-line determination and elimination. 
Limitations prior to that time had been on the types and shapes of 
objects considered. With point visibility schemes applied to complex 
objects, calculation times had increased to unacceptable levels for many 
applications. Appel introduced his idea of "quantitative invisibility" 
to reconcile surface and point visibility approaches. 

Freeman and Loutrel [Ref. 4] were concerned with the two-dimensional 
version of the general hidden-line problem. This limitation was imposed 
in order to work with arbitrary, irregular polygonal figures. 

Comba [Ref. 2] addressed the problem of detecting intersections of 
three-dimensional objects mentioned earlier by Johnson [Ref. 8]. 
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Included was a discussion of the overall hidden-line problem as well as 
several interrelated sub-problems. Prior efforts by other investigators 
were mentioned. 

A recent contribution by Matsushita [Ref. 10] used initial hidden- 
line computations for subsequent calculations involving rotation about a 
single axis. Reduction of processing time requirements was the result. 

As seen in earlier research, computation time constraints were proving 
to be a major motivating force for work undertaken in pursuit of problem 
Solutions to real-time applications. 

The reader is directed to a comprehensive review of perspective view 
and hidden-line elimination considerations by Desens [Ref. 3]. Included 
in the thesis was an extensive bibliography. 

A pair of articles dealing with hardware aspects in relation to the 
problem areas at hand should be mentioned. 

Hagan and Treiber [Ref. 6] provided an interesting approach to 
generalized linear transformation, spatial coordinate transformation, 
generation of visual displays for three-dimensional objects in perspect- 
ive, generation of functions of two or more variables, and convolution 
and correlation operations. The authors described a number of hybrid 
analog-digital techniques using multiple analog and hybrid elements 
in parallel. 

A hybrid, parallel array, composed of digi tal -to-analog converters 
(DAC's) and hybrid multiplying elements (HME's), were used for rotation, 
translation, and perspective representation of an object. This spatial 
coordinate transformation was performed on each point defining an object. 
The authors stated that using an analog vector generator which develops 
a line segment in 30 microseconds, it was possible to display objects 
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consisting of about 1000 line segments at a frame repetition rate of 
30 frames per second. 

It was necessary to perform nine multiplications, 12 additions, and 
two divisions for each operand fetch of the values representing a point. 

The parallel nature of array elements reduced the number of arithmetic 
operations per operand fetch to slightly more than seven. These developed 
techniques had proven useful for a wide range of real-time, signal 
processing applications. 

Hagan, Nixon, and Schaefer [Ref. 7] presented background motivation 
for design of the Adage line of graphics terminals. Rotation, depth 
cueing, and windowing techniques were specifically discussed in terms of 
hardware implementation. Basic software operators were also mentioned. 
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II. EXPERIMENTAL PROCEDURE AND RESULTS 



Experimentation commenced with programming to implement a graphical 
display, following reference research into previous work done on the 
problem to be investigated here. 

A. GRAPHICAL DISPLAY OF THREE-DIMENSIONAL OBJECTS 

Programming and test runs to allow display of three-dimensional 
objects were attempted first. The concepts discussed and techniques 
developed by Desens [Ref. 3] formed the basis of a program. One of the 
basic goals of the program was to transform the rectangular coordinates 
of an object into two-dimensional perspective for CRT display. 

A single 4x4 matrix, [T], performs rotation, scaling, translation, 
and perspective transformations. Partitioning of [T] separates those 
four operations. The upper left is a 3 x 3 rotation matrix. The first 
three elements of the bottom row hold the translation terms. The first 
three elements of the right column contain inverses of the three viewing 
distances necessary for perspective transformation. The bottom right 
element is a single scale factor. 

The matrix operates on a homogeneous coordinate vector. Original, 
orthogonal coordinates (x,y,z) become the homogeneous coordinates 
^^h’^h’^h’^h^ where x^ = Sj^x, y^ = Sj^y, and = s^^z. The added, fourth 
coordinate, S|^, is considered a variable scale factor. 

Transformation of a point is made by post-multiplying a homogeneous 
coordinate vector by the complete transformation matrix, [xj^,y|^,Z|^,s^] 
[T] = ,y|^ ,Z|!| ,S|!|] . To obtain the final display coordinates for the 
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X-Y plane, division of the transformed homogeneous coordinates by the 
S|^ coordinate is made, that is X = X|^/S|I| and Y = 

The program included a two-stage procedure for removal of hidden- 
lines. First, all lines hidden by an object's own volume were removed. 
Next, each segment not yet eliminated was examined over it's length and 
that portion remaining in view was displayed. 

Two test objects were used to supply data points for experimental 
computations. An L-shaped block defined by 12 points was first displayed. 
The second test object, an aircraft carrier, consisted of 62 points. 

Each data point was specified in rectangular coordinates. 

The results from tests were successful in terms of the resulting 
graphics display output. However, the preliminary results were not 
satisfactory from two standpoints. While the perspective view trans- 
formations provided the proper depth cueing, the hidden-line elimination 
segments of the program used were not functioning currectly. Addition- 
ally, processing time requirements were excessive for a simple and a 
moderately complex object, i.e., the two test objects. 

Confirmation of Desens' conclusion [Ref. 3], regarding the relation- 
ship between processing time requirements and picture complexity, was 
made on the basis of the test runs. Picture complexity was defined by 
the number of plane-faced surfaces included in an object and thus the 
number of lines defining those surfaces. The conclusion validated by 
the test runs was that time required for connection of points increases 
rapidly as picture complexity increases, approximately in proportion 
to the square of the number of lines. 
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B. HIDDEN-LINE ELIMINATION 

A different approach to the hidden-line removal problem was taken 
next. A simple, fast algorithm with emphasis on hidden-line elimination 
was sought. Attention was directed to the' hidden-line plotting program 
developed by Williamson [Ref. 22]. 

1 . The Basic Program 

Subroutine HIDE, Appendix C, yielded a two-dimensional representa- 
tion of a surface or object. Curves or segments of curves were plotted 
where not hidden by segments previously plotted. Any curve segment lying 
above a visual maximum function was determined to be visible. 

Options included in the program were translation of data points 
to simulate stepping in the depth dimension, plotting of a border, axes, 
and a title, and plotting of the visible portion of the underside of a 
surface. In the latter case, segments of curves were presumed to be 
hidden where they lay above those segments previously plotted. Provisions 
for perspective view or rotation transformations were not included in 
the program. The appropriate transformations could be applied prior to 
a call to subroutine HIDE. 

The first call to HIDE was for initialization and output of the 
curve farthest in the foreground. Subsequent calls to subroutine HIDE 
performed hidden-line determinations and output of visible curve segments 
farther in the background. 

The main program computes the x and y coordinates of the curves 
(function) to be displayed. Additionally, the calling program sets 
initial values for control of display options. 



15 



2. Extensions 



The basic program was to be extended to encompass graphical 
display output and as many interactive options as possible within the 
limits imposed by time constraints, 
a. Graphical Display 

Adaptation of the program to graphics terminal display was 
the first step. Equipment to be used consisted of an Adage Graphics 
Terminal (AGT-10) and an SDS-9300 computer, both located in the Naval 
Postgraduate School Computer Laboratory. 

The SDS-9300 is a medium sized, general purpose, digital 
computer system. The sophisticated central processing unit can handle 
a large number of operations. Main memory consists of 32K words of core. 
A magnetic drum is used for secondary storage. The system is controlled 
at the operator's console and at a teletypewriter. Input is via card 
reader, teletypewriter, paper tape reader, or magnetic tape. Output is 
sent to a line printer, paper tape punch, teletypewriter, or magnetic 
tape. Communication links are available to two Adage AGT-lO's. 

The AGT-lO's are small, general purpose, digital computers 
with each having 4K words of main memory plus a magnetic disk unit for 
secondary storage. A graphics display console is linked to each AGT-10 
for text and graphics presentation. Inputs to an AGT-10 are a tele- 
typewriter, paper tape reader, lightpen, joystick, function switches, 
and variable control dials. Output appears on a CRT, paper tape punch, 
or teletypewriter. 

Display capability was implemented through the use of systems 
graphics routines DGINIT [Ref. 16], IHEAD [Ref. 12], IPACK [Ref. 13], 
and GRAPHO [Ref. 18]. The purpose of FORTRAN callable subroutine DGINIT 
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is to initialize the AGT-10 graphics display subsystem. DGINIT must be 
called before any graphics communication may occur between the SDS-9300 
computer and an AGT-10. Function IHEAD constructs the one-word header 
of a graphic data block to be displayed by a subsequent call to GRAPHO. 
One of the two parameters passed to IHEAD determines whether all the 
vectors drawn are to be solid lines or dashed lines. The other parameter 
specifies display intensity. FORTRAN callable function IPACKwill pack 
the X and y values and a move-draw indicator for a coordinate point into 
one word of a graphic data block. Subroutine GRAPHO transmits a graphic 
data block from the SDS-9300 computer to an AGT-10 for display. 

Scaling of data points was necessary for reconciliation of 
differences in the real number values used by the graphics terminal and 
the actual size or scale of the displayed object. 

In order to maintain a reference point from which to work, 
data points were maintained in their original, unsealed format. A 
parallel data structure, array FIGURE, was created to provide the vehicle 
for storage of the scaled, packed data words. 

The communication link between the SDS-9300 computer and an 
AGT-10 graphics terminal was established through AGT-10 systems program 
GATED. Before a user's program is executed on the SDS-9300 computer, 
GATED must be loaded and executed on an AGT-10. The purpose of GATED is 
to allow display (and optional editing) of graphics and/or text data 
blocks received from the SDS-9300. GATED is further described in more 
detail in Ref. 17. 

The initial scheme for determining which data points were to 
be displayed consisted of witholding display of visible portions of a 
curve as they were computed until the visual maximum function had been 
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completely updated for the current curve, i.e., current call to sub- 
routine HIDE. At that point the complete visual maximum function was 
transmitted to an AGT-10 for display. 

Display in this manner, upon succeeding calls to subroutine 
HIDE, was unacceptable. Where hidden lines occured, portions of the 
visual maximum function were displayed more than once. This situation 
caused an accompanying increase in the intensity of certain curve seg- 
ments. The gain in intensity was directly proportional to the number 
of times the same curve segment was displayed as a portion of the visual 
maximum function. 

Changes were incorporated in the program to allow graphical 
display of only those portions of a curve lying above the visual 
maximum function. In this manner each visible curve segment was displayed 
only once, thus solving the problem of non-uniform intensities. However, 
a problem of storage requirements developed. 

The interface systems program, GATED [Ref. 17], limited the 
number of graphics blocks allocated for display to 51. That was found 
too restrictive for the number of visible curve segments encountered in 
most aspects of the test function display. Anticipated requirements for 
future applications also exceeded the existing constraints imposed by 
GATED. 

Consequent modification of GATED included allowance for more 
graphics blocks. The increase in storage required by the greater allow- 
ance of graphics blocks was more than offset by a coincident decrease in 
the number of text blocks. 

An alternate solution to the problem of an insufficient 
number of graphics blocks for display of the entire test function involved 
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programming changes within subroutine HIDE. Instead of displaying each 
visible curve segment after being computed, display was witheld until 
all visible portions of a single curve were determined. The data points, 
representing the visible curve segments, were packed into consecutive 
locations in array FIGURE. Discontinuities between endpoints of segments 
were accounted for with a "move" indicator in function IPACK [Ref. 13]. 

A "move" rather than a "draw" was indicated for the vector from the last 
point of a visible segment to the first point in the next visible seg- 
ment. Once all the segments were computed and packed, array FIGURE was 
displayed as a single graphics block. There then existed a one-to-one 
correspondence between the number of curves and the number of graphics 
blocks required. Visual results showed all the visible portions of a 
curve at the same instant. The preceding method had displayed each 
segment separately. 

Consideration was given to display of the entire series of 
curves, representing the test function, at one time. This method would 
have involved, prior to display, storage of all visible points for the 
entire structure. The scheme used required temporary storage for display 
of each curve separately. While storage of this magnitude was avail- 
able in the SDS-9300 computer, utilization was considered extravagant. 
Local variable storage requirements, in the case of the test function, 
would have increased to a total of approximately 12,000 words versus 
the 3,000 words required to include storage of one curve at a time. The 
increased storage would have been manifested in a highly expanded array 
FIGURE. 
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b. Test Function 



The previously mentioned constraint imposed by subroutine 
GATED on the number of graphic data blocks also affected the selection 
of a suitable test function. The curves were consequently limited to 
100 points each. This limitation enabled display of the test function 
with no discernible degradation in the smoothness of the curves. The 
number of curves chosen to represent a surface was 30. When this number 
of curves was combined with the appropriate translation values to simu- 
late stepping in the depth dimension, a reasonable representation of a 
curved surface resulted, as in Figure 2. The number of data points, 
required by both the top and underside of the test function, caused 
a core overflow condition in the AGT-10. Additionally, a second call 
to HIDE was deemed to be unnecessary for the possible display features 
desired by the problem statement. For a series of two-dimensional 
curves offset for simulated depth stepping, i.e., the signal processing 
example (Figure 1), it was expected that coordinate points would be 
specified for positive half-space only. 

Deletion of the second call to subroutine HIDE resulted in 
a significant decrease in processing time, since for even the hidden 
portions of the underside of the test function the majority of calculations 
in subroutine HIDE were made. 

c. Interactive Extensions 

Man-Machine interaction is a necessary and sufficient condi- 
tion to justify the use of a graphics terminal. On-line modification 
of a display and program control flow options exercised by an operator 
sitting at a graphics terminal are highly desirable in many applications. 
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Figure 2. Test Function 

(1) Variable Control Dials . The control dials were chosen 
as the method of input since they provided a range of values, not just 
simple switches. Also, they would maintain their settings until changed 
by the operator. 

At the time experimentation was being conducted, no 
method existed for making available to the SDS-9300 computer, and thus 
to the FORTRAN program, the values of the variable control dials. While 
not a primary consideration of investigative work, control dial implemen- 
tation resulted in a significant part of resource expenditures, time 
and effort. 

The development of a utility program to use the control 
dial settings was accomplished. Appendix D contains a discussion of 
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implementation details and programming modifications. Appendix B 
includes instructions for on-line display modification and option 
selection using the variable control dials. 

(2) Translation . A simple static translation operation was 
applied to the data points after scaling required for the graphics 
terminal display and prior to being sent across the communications inter- 
face to an AGT-10. This operation allowed an operator freedom in place- 
ment of a displayed structure on a CRT. Operator control of placement 
was input via variable control dial settings. 

(3) Rotation . Rotation of the display structure about two 
axes was simulated. An apparent rotation of the surface represented 
about the y axis resulted from translating all points of succeeding' 
curves in a negative x direction. In a similar manner, translation of 
all points of succeeding curves in a positive y direction resulted in 
simulated rotation of the surface about the x axis into positive half- 
space. Referring to Figure 3, the receding curves have been rotated 
(translated) to the left and up from the orientation in Figure 2. 

Translation operations to simulate rotation were applied 
to arrays' XH and H in conjunction with other modifications prior to 
graphical display. XH and H are working arrays initialized to the current 
visual maximum function in subroutine HIDE. 

The apparent rotation effects produced by translations 
v;ere acceptable for small to medium sized angles, i.e,, up to approxi- 
mately 40 degrees about the two axes. For larger angles, distortion of 
the displayed surface was seen. Part of the distortion could be explained 
by the fact that the curved surface was represented by a series of curves 
displayed close together and as the apparent rotation angle was increased. 
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the curves became more disjoint thus destroying the illusion of a 
continuous surface. 




Figure 3. Rotation of Test Function. 

Applications specifying a two-dimensional series of 
curves would not suffer from the display limitations imposed by this 
particular scheme as would three-dimensional object representation. For 
a series of curves, variable offset in the apparent viewing angle would 
be provided for possible qualitative trend analysis. The operator again 
had control over the variance of rotation (translation) by adjusting 
control dial settings. 

(4) Window Operation . To enhance the capabilities of the 
developing program and to enable an operator to improve any qualitative 
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analysis of displays, a windowing operation was added. The window was 
envisioned as a planar slice through the displayed surface or family 
of curves. All points of the plane would have the same y value in 
three-space, i.e., parallel to a plane defined by the x and 2 axes 
(Figure 4) . 




Figure 4. Windowing Operation. 

In terms of programming, synthesis of the window was 
accomplished by merging a straight line, representing the window in 
relation to one curve from a series of curves, with the previous visual 
maximum function. A new visual maximum function was thus formed and was 
used for computing visible curve segments. The window acted as a lower 
threshold level , below which no curve segments were displayed. The 
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window "line" is translated (updated) for successive curves to provide 
the illusion of a planar slice in the same orientation as the display 
structure. 

Modification was made to the program to include display 
of the first curve after hidden-line calculations were made. The previous 
concept had been that the first curve defined the initial visual maximum 
function for use with succeeding curves. With the introduction of a 
window operation, that concept was not valid. In other words, a portion 
of the first curve may also be hidden by the window (Figure 4). Operator 
control of the window was implemented with variable control dials. 

(5) Cursor . Provision for some method of quantitative 
analysis was desired. Addition of a cursor was the next goal. As 
envisioned, the operator would be able to select the cursor for any of 
the series of curves. The cursor should be adjustable along the ordinate 
and of varying, controlable length to specify points on a curve for 
measurement. Also, placement of the cursor at points other than on a 
curve might be useful for determining relative points on different curves. 
The operator should have the option of suppressing, at any time, a cursor 
mode of operation, since the continuation of successive curves being 
displayed would be held pending the end of cursor analysis. The accuracy 
of data point location by cursor placement was not expected to be of a 
high degree due primarily to discrimination limitations. Operator 
matching of cursor and curve points was a crude, viusal process. 

The cursor, as implemented for display, consisted of a 
single vector defined by two endpoints. The first, or reference, end- 
point had a y value corresponding to a value of zero for the particular 
curve being analyzed. The second, or variable, endpoint had an operator 
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controlled y value. The x values for Doth endpoints were equal and both 
were determined by the operator. 

In Figure 5, the reference endpoint is at the top and 
the variable endpoint has been adjusted to lie on the last curve 
displayed. The variable endpoint position represents a negative y 
'value, lying below the reference endpoint. 




Figure 5. Cursor with Text. 

The cursor could be invoked by an operator using a 
particular setting for one variable control dial. Once the cursor was 
displayed, a different dial controlled the ordinate values, and there- 
fore the position on the CRT of the cursor. Another dial determined 
the abscissa value of the variable cursor endpoint, and therefore the 
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length of the displayed cursor. Still another dial value was tested in 
the program to determine whether to continue in the cursor mode for the 
current curve or to resume display of succeeding curves. 

Display of the next curve caused the cursor to be blanked 
from the CRT since the next curve used the same graphics data block number 
as that used by the cursor. Output of the x and y values of the variable 
cursor endpoint were written’ on a line printer. The option to print the 
values was provided to the operator through another variable control dial. 
A call to systems subroutine DELAY was necessary to prevent the same 
values from being printed several times while an operator was modifying 
the dial settings. 

While the basic implementation of the cursor seemed 
acceptable in terms of analysis capabilities, the remote output of vari- 
able cursor endpoint values was not adequate. Additionally, control 
dial manipulations were awkward. 

Modifications to subroutine CURSOR were undertaken to 
eliminate the two deficiencies mentioned. Textual display of variable 
cursor endpoint values, appropriately scaled for the curve, would provide 
instantaneous feedback to the operator. The interactive nature of the 
analysis capability of the program would thus be increased. 

Addition of a line number, displayed to aid the operator 
by providing a reference, was made. Display of the x and y variable 
cursor endpoint values were included in the same programming loop used 
to sense the control dial settings. By being in the same loop, only 
the value portions of the text were changeable once subroutine CURSOR 
was entered. Figure 5 is an example of cursor operation with accompanying 
text. 
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System subroutines used for text display were DTINIT 
and TEXTO [Refs. 11 and 15]. DTINIT initializes the AGT-10 text display 
subsystem and must be called before any text communications may be 
initiated between the SDS-9300 computer and an AGT-10. TEXTO allows 
output of a full (or partial) line of text on an AGT-10. Prior to a 
call to TEXTO, text characters are packed into a binary coded decimal, 
text array with four characters per word. The text array is then passed 
as one of several parameters to subroutine TEXTO. Formats for encoding 
of text characters are the standard FORTRAN I/O formats. 

Return from the cursor mode of operation was now 
controlled by the adjustment of only one dial. Operator intervention, 
in terms of control dial manipulations, was significantly reduced. 
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III. CONCLUSIONS 



Following is a discussion of summarized results and program 
extensions with potential applications in mind. 

A. SUMMARY OF RESULTS 

Graphical display implementation of a hidden-line removal algorithm 
has been accomplished. The program produces a two-dimensional represen- 
tation (CRT display) of a three-dimensional surface or series of curves 
by displaying segments of a succession of curves where the segments are 
not hidden by any of the curves previously displayed. 

Capabilities of the basic program were extended to enhance the design 
and analysis potential for applications use. Interactive attributes of 
the existing program include translation of the display structure for 
positioning, translation to simulate a limited rotational capability for 
presentation of different viewing angles, a windowing operation for 
designating threshold values of displayed curve segments, and a cursor 
mode of operation for analysis. 

The operator is able to affect his desired changes in the display by 
input of variable control dial settings. An extended version of GATED, 
the systems communication interface routine, was implemented as P0G3 to 
allow the operator this method of input. 

B. DEVELOPMENT AND APPLICATIONS 

Any extensions to the program, manifested as user options, will 
primarily be a function of specific applications. Reinclusion of a 
second call to subroutine HIDE will allow display of the unhidden part 
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of the underside of a curved surface. Changes to subroutine HIDE, as 
listed in Appendix C, will be necessary for consideration of the sign 
of the abscissa values displayed. The reader is cautioned to consider 
data point storage limitations and processing time increases caused by 
a second call to HIDE. 

The option of displaying a border, axes, or title may also be desired. 
Scaled axes would provide a reference framework for estimating curve 
point values and orientation of the display as viewed. 

Additional rotation capability is desirable. Rotation into quadrants, 
other than those presently allowed, will extend the present limitations 
imposed on viewing angles. The problem of figure distortion for large 
angles will, however, remain. • 

Calculations to result in perspective view display would significantly 
increase the class of possible applications. An example of the types of 
calculations necessary can be found in the transformation schemes used 
by Desens [Ref. 3] and Johnson [Ref. 8]. 

No consideration was previously given to text or graphics editing. 

These features are available with certain display devices. Editing of 
displayed structures may well be necessary for certain design applications. 

As an alternative to design using editing features, it might be 
advantageous to allow the program user the ability to change parameters 
of a function, thereby modifying or creating his own figures. This 
capability could be implemented with additional control dials or similar 
input devices. Associated with parameter variation might be a "menu" 
selection option. The user would specify a basic function or object, 
for example with a lightpen, for later design modification through 
parameter variation. . 
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The potential applications of hidden-line removal in the areas of 
design and analysis are many as amply documented in several of the 
references contained in the bibliography section. A simple, fast 
algorithm, as adapted here, is suited to a wide range of problems that 
call for display of real-time solutions. Interactive, man-machine 
communication is another essential ingredient of real-time applications. 
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APPENDIX A 



FLOWCHARTS 




FLOWCHART 1 
MAIN PROGRAM 
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APPENDIX B 



USING THE PROGRAM 

A. EXPLANATION OF VARIABLES 
I . Main Program 



X 


The ordinate array for one curve. 


Y 


The abscissa array for one curve. 


Y1 


A working array used for intermediate values in 
computing Y values. 


XG 


The ordinate array of the visual maximum function. 


G 


The abscissa array of the visual maximum function 


XH 


A working, ordinate array used in conjunction with XG 


H 


A working, abscissa array used with G. 


IGDIR 


The first word address of the graphics directory used 
as a parameter in calling DGINIT. 


DIAL 


The array used to hold unsealed values of the 
variable control dials. 


MAXDIM 


The dimension of arrays XG, G, XH, and H. MAXDIM 
is used for overflow condition testing of those 
arrays. 


XMIN 


The lower bound for X values. 


DELTAX 


The ordinate data increment. 


YMIN 


The lower bound for Y values. 


DELTAY 


The abscissa data increment. 


XLNTH 


The horizontal scaling factor used in computing 
translation. 
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YLNTH 


The vertical scaling factor similarly used for 
translation. 


XX 


A constant used in computing Y values of the test 
function. 


IDEV 


Specifies which graphics terminal is to be used 
for display, either 1 or 2. 


STEP 


Has two meanings in main program. First, STEP is 
used as the ordinate value increment. Secondly, it 
is used as an incremental value in determining 
translation of succeeding curves. 


TX 


The scaled ordinate translation value for 
positioning of display 


TY 


The scaled abscissa value for positioning 
translation. 


YWINDO 


. The scaled value of the window operator. 


NG 


The number of points in the current visual maximum 
function after each return from subroutine HIDE. 


N 


The number of points defining a curve to be 
connected by vectors for display. 


NFNS 


The number of curves defining a surface to be 
displayed. 


NUMFIG 


An index value used to keep track of sequential 
graphic data block numbers for calls to GRAPHO. 



2. Subroutine HIDE 



HDIAL 


The array for local values of the variable control 
dials. 
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EPS! . 



The relative abscissa increment used to simulate 



F(XX,XI 

DXIN 

DYIN 

FNSMl 

ABS(a) 

INDEX! 

YZ 

AMAX(a, 

EPS 



discontinuities in the visual maximum function. 

EPSl should be one or two orders of magnitude larger 
than the smallest recognizable relative difference 
in single-precision, floating-point arithmetic. 

,YI ,XIP1 jYIPl ) A statement function used to compute 
the ordinate value on the line joining the point 
(XI, YI) and the point (XIPl , YIPl) corresponding 

to the abscissa XX. - 

The incremental ordinate translation value to 
simulate stepping in the depth dimension. 

The abscissa translation value similarly for 
simulated depth stepping. 

The number of curves to be translated. 

An intrinsic system function to calculate the 
absolute value of real number a. 

An index used for initializing the visual maximum 
function. Later used as an index during X-table 
lookups. 

A temporary variable used to compute the initial 
visual maximum function. 

)) An intrinsic system function used to determine 
the largest of two real numbers a and b. 

The incremental value added and subtracted from 
• curve endpoint ordinate values. 
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DXKK 


The X translation value to simulate depth stepping. 
Each curve is translated relative to the first curve 
displayed. 


DYKK 


Similarly, the y translation value for depth 
stepping. 


RELINC 


The slope of the x,y translation vector. 


INDEXG 


An index used for lookup of the XG table. 


LAST 


A flag for signaling that the end of X-table has 
been reached during lookup. 


FIGURE 


The array to hold scaled, packed data points for 
display. 


IWHICH 


A flag used in interval evaluation for hidden-line 
determinations. 


SLOPE 


The slope of the vector joining endpoints of the 
interval in the region where curve crosses visual 
maximum function. 


NGRAPH 


The number of points defining each visible curve 
segment. 


NWD 


The number of data words defining all visible 
portions of a single curve. 



B. OPERATING INSTRUCTIONS 
1 . Graphics Terminal 

The following should be entered on the terminal teletypewriter 
to execute P0G3: 

START("P0G3",pvv)! 

GATED! 



where pvv is the pack and volume number location of P0G3. 



2. SDS-9300 Computer 

The SENSE 2 button on the console should be depressed prior to 
reading in the card deck. This assures proper interpretation of the 
patch deck at the beginning of the program. 

Immediately after the program starts execution, the user will 
notice an input light on the teletypewriter adjacent to the computer 
console. The proper operator response is: 

IDEV=1* (carriage return) 

This specifies AGT-10 number 1 for display. At present the variable 
control dials are available for use with unit number 1 only. 

3. Variable Control Dials 
a. Initial Settings 

The following are suggested initial dial settings for use 
with the test function as listed in Appendix C: 



DIAL VALUE 
A and D 0.33 
B and E -0.33 



RESULT 

Determines an initial orientation. 
Positions the displayed structure 
approximately about the center of 



the CRT. 



C 



F 



0.00 to Allows continuance of succeeding curve 

0.10 display. 

-0.99 Sets window to maximum negative 

value to suppress it's effects on the 
initial test function display. 



b. View Orientation 

To simulate rotation about the y axis, adjust DIAL A, 
positive values. 
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To simulate rotation about the x axis, adjust DIAL D, 
positive values. 

c. Display Positioning 

To modify horizontal (x) position, adjust DIAL B. To 
modify vertical (y) position, adjust DIAL E. 

• d. Window 

Increasing the value of the window can be accomplished by 
increasing the setting of DIAL F. 
e. Cursor 

To display the cursor and accompanying text for a particular 
curve, set DIAL C greater than 0.10 immediately prior to display of the 
curve in sequence. 

To modify x coordinates of cursor, adjust DIAL A. 

To modify y coordinate of cursor variable endpoint, adjust 

DIAL D. 

To exit from cursor mode and resume display of succeeding 
curves, set DIAL C less than or equal to 0.10. 
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APPENDIX C 



COMPUTER PROGRAM 
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CALL Lf SKUP(ZZ/X( rOEXT )><1 ) 
CALL LC0'<'JP(7Z>XG( I'JDEXG)i<2) 
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YF = H(Y)/5. + TY 
FIGUI^E(IZ) = IPACK(XF/YFi 1) 



o 



iD 



O 



II II H- 



rx! 



X 
II LiJ 

































o 






. — - 


























d: 










o 






4—4 


























LjJ 










IN 
































t — i 
















+ 


























•d 










CD 
































CO 


CO 




4— » 




h- 






CO 


























4 


►—4 














CO 


























Li. 


Ll 




1 




3 






.-4 


























✓ 


3T 








CO 
































Z) 


3 




CO 










CO 




























X 














X 


























-d 










r 


































ccr 




II 




• 4 
































o 


Ll! 






fO 


Q 






CO 




























*~4 




CO cn 


X 






o 
































-p' 




<r 






►—4 


























■ — 


— 






cn 


>■ 






— ■ 


























^ — 1 


•m — 1 






h- 


• 


CO 




CO 














O 














o 




■ — 




L- 






% 






^ — i 








O 












LiJ 


« — 1 




4— ♦ 


cz 


O 


Ll 










nO 














Oi 






(T 








(0 


• 






CO 














■rH 






vO 






3 






• 




CO 


+ 




CO 






(Tj 




















CO 


1— 






. — . 


CO 












K- 








cn 






<r:> 






•— t 


3 


4 — 5 


o 




•—4 


. — . 




— 














h- 












u. 


X 




X 


4— e 




3 




CO 




CO 


O 






















H- 


+ 






1 


3 




X 






CO 








(D 






cn 






> 


30 




CO 


■ — 






















CO 






CO o 




LU 


(D 


CO 


X 


X 


CO 


X 








T— i 


















r-i 




O 




K—4 


• 


» 




X 




4—1 




CL 




















CO 


*— 




Ll 


Uj 


Uj 








X 


+ 


CO 


• 








^ — 1 






T-H 


INI 




’ — ■ 


O 


y_ 


3 


3 


-L 


II 




— 




CO 


C.^ 








T' 






• 


•— 


II 


CD 


* 


“O 


m 


• 






•4“< 


X 


CO 




Ul 








• 


00 


1 


LiJ 






X 








— 


m 








CO 




• 








u> o 




_J 


II 




X 








c? 




^ — 1 


+ 


u 


♦—4 


II 


H- 








_J o 


N1 


• 




■n 


< 


• 


II 


















co 


ai 


OJ 


a 


• 




t— 


o 


X) 






d: 






■ — 




+ 


3 


II 


II 


3 


c 


X 


u 


INI 


1— 










LJ 




'uJ 


CO 


CO 


LO 


3 




3 








-J 










Cl 


II 


~y^ 


00 


CC 




•~4 


►-4 


X 


X 


3 


3 




—V 


4—1 




' — ' 


II 


II 


II 


' — • 


h- 






o 


10 


3 




u. 




— 


— 


3 


II 


3 


CL 


3 




















to 


_j 




z. 








■ — 




3 


cn 




Lu 




,r— 1 


n — ^ 


Ll 


(T» 


■< 


Ll 


CD 


►-4 


c 


Ll 


3 


Ll 


Ll 


Ll 


”0 


3 




CO 


O 




X 


u. 


M 




LZ 




►-4 


O U 


LJ 


►-4 




•— * 


►-4 




X 


3 


X 


►-4 


O 



— • C'T “3 Oi 
X — 

+ C -f 

II 

0 II ^ 

^ ~y 



; X X 



^ o 
o o 
o o 



o 

vC 



sD 



m 

vO 



OJ 

vO 



(n 

cn cn 



m 



55 
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YWIN30 = Yl'. I\DB + DYIM 
DO 13 1 = l.JJ 
G( I ) = AYAX ( YWI\'D6/H( I ) ) 

XG( I ) = Xr)( I ) 



o 

•— I 

u 

5 : 

3 

•s 

- >- 

OJ o 

o 

u 

> M 

X 

o 

- cr — 

•V o X 

a: c/) »~ 

LiJ C:T. -s 
—I 3 X 
CJ H- 



C/) 














3 


> 










cu 












0 


3 


lU 










UJ V 


m 






V 


_J 




< 


0 










V 


in 






V 


< 




CJ 












+ X 








X 


















0 


CD 






0 


0 


X 
















X 








X 


3 


X 1 












“:> + 






X 


1 




CL 


• 












“> 


<D 




X 




q: 


X 


0 












^ 


CD 




X 


— 


bJ 


3 


• 






x: 






CD — * 






0 




►— t 


Cj 


X 






►— < 






X >: + 


— 


X — 1 




— 






CD 






Q 






X 


0 


7^ 


-1' 


X 


> 




» 






X 






II ”3 


• 


•s 




X 


Ui 


0 








< 






II “) 


X 


X — { 




~p^ 


0 


• 








> 






' 


_J 




» — « 


CD 


*— t 


y 


— 






1 






^ ^ <C. 


• 


It 


— 


• — 




^-. 


3 












X— • 


CO 




X 


00 


0 


% 


C 






II 


vO 




+ II 


71 


•— 








cr 


t — 1 








vO 




+ 


U 




II 


II 


> 


u» 0 




T' 


5: 










OvJ 










X 




cr 




CD 




”) ^ C3 


— 


00 






_J 








3 


0 X 














— t 








X 


X 




0 


0 -- -- 


Lu 


<D 


■ — 


— 


<t 


u_ 


Li- 




UJ 


<r 


Co 




X a CD 


.—• 


0 


X 


X 


0 






X — f 


cr 


5: 


CD 


UJ 




vO 






OJ 


ro 










0 








vO 






00 


LO 










0 







56 



subroutine L09KUP(X/XTBL> J) 
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FMC03E (8/60/ I T/T2 ) X 
60 rt)^^.AT(F8.8) 

CALL TfXTGC I DEV/ I TXT?/ ?/ 38 / 5/ 1 / 3/ I ER ) 
FNCUDE ( 8 / 7"'/ ] TXT? ) Y 
70 re<8 yat ( r 8.? ) 



CALL TLXT9( I DEV/ I TXT?/ 2/ 40/ 5/ 1 / 3/ I ER ) 
IF (CDIAL(3) .GT.0.1 ) 33 T3 99 
CALL 0TInIT( IDFV/ ITOIR/ 7/ lER) 

RETURN 

END 
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LOA =0 

B«U VCDB 

VCODEr^ LDA =1 

VCCL’-^ VCD 

VCDSWe P 7 .E 0 



PZt 

PZE 

PZE 

END 
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APPENDIX D 



CONTROL DIAL IMPLEMENTATION 

An initial change was made to GATED [Ref. 17], the AGT-10 program 
to display and edit graphics and text data blocks. The alteration in- 
cluded a call to AGT-10 system subroutine TRVCD [Ref. 19], which obtains 
the digitized values of the output voltages of the variable control dials. 
The sampling and storage of the six-dial values requires 750 micro-seconds. 
The sub-program, which includes the call to TRVCD, was titled POGG. It 
thus became a required subprogram of P0G3, the revised version of GATED. 

In addition to sampling the control dial settings, POGG transmitted 
the six values to the SDS-9300 computer and stored them in upper-core 
locations. POGG was invoked once every cycle through P0G3. This meant 
that current dial values were available every 0.025 seconds. 

Once the dial values were resident in the 9300 computer, they were 
transferred to locations designated by a FORTRAN array name, e.g., 

DIAL. At this point scaling of the dial setting was completed and the 
resultant values passed to subroutine HIDE for display purposes. 

An intermittent difficulty occurred with control dial testing. The 
problem manifested itself as sudden and apparently random blanking of 
CRT display. Eventual diagnosis was interference between sending control 
dial information, from the AGT to the 9300 computer, and coincident trans- 
mission from the computer to the graphics terminal of a graphics display 
request, GRAPHO [Ref. 18]. 

Solution of the interference problem involved the use of a meta- 
symbol subroutine, VCD. The purpose of VCD was to limit the sending of 
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control dial values to the SDS-9300 computer. Values were then sent only 
after VCD had been called. This interrupt condition insured that no 
conflicts would occur. 

The inclusion of the metasymbol subroutine VCD forced further modi- 
ficiations of P0G3. A patch deck was necessary at the beginning of the 
program. The patch deck used is included in the program listing. 

Appendix C. A test for code 35, the designation given to a request for 
control dial sensing in P0G3, was included. Additionally, changes to 
instructions in AXINT, a required sub-program of P0G3, were made to pre- 
clude a jump to the magnetic tape operation. This operation had the 
only other 30 series request code. A SETGT instruction was substituted 
for several consecutive operators in AXINT to allow for code 35. 

P0G3 is currently available on the user file pack for AGT-1 in the 
Naval Postgraduate School Computer Laboratory. Instructions for loading 
and executing P0G3 are included in Appendix B. Following is a summary 
of programming changes to GATED. 

1. Title changed to P0G3. 

2. After third instruction under label SA3 and immediatly after: 

JPAN SB 

the following code was added: 



MDAR'H 


$GTWD1 


ARRS 


3 


MDAR'F'A 


77 


MDXO'F 


35 


JPLS 


.+5 


JPSR 


$P0GG 


JPSR 


FSWI 


. 


VBLK 


JUMP 


SB 



The code sequence resumes as in GATED with 
MDAR'H $GTWD1 
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3. After label 6BLK: 200 !H, add 

VBLK: 350 IH 

4. File program POGG required by P0G3 is as follov/s 



EXPUNGE 






TITLE 


POGG 




ENTRY 


POGG 




POGG: 


JUMP 






JPSR 

ARXO'F 


$TRVCD 




ARMD 


CNT 




MDAR'F 


$TVCDA 




MDAE'L; 


-1 


FADD: 


0 






MDAR'X'I 


FADD 




ARRS 

NOOP 


6 




ARMD 'I 


FADD 




MDAR'X 


CNT 




MDXO'L; 


6 




JPLS 


FADD+1 




JPSR 


$ROWFW 

-0 






POGGl 






$TVCDA 

6 




MDIR 


POGG 


POGGl : 


77750 




CNT: 

TERMINATE 


0 
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